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Abstract. Dust configurations play an important role in astrophysics and are the simplest models 
for rotating bodies. The physical properties of the general-relativistic global solution for the rigidly 
rotating disk of dust, which has been found recently as the solution of a boundary value problem, 
are discussed. 



1 Introduction 



Dust, as the simplest phenomenological material, is a good model for astrophysical and 
cosmological studies. As a source of gravitational fields, dust may be interpreted, in a hy- 
drodynamical language, as a many-particle system, whose particles (mass-elements) interact 
via gravitational forces alone. While the cosmological relevance of that model has been de- 
monstrated very early p], investigations of isolated dust configurations are rather rare [2j. 

This paper is meant to discuss dust as a model for rotating bodies and to focus attention 
on the relativistic and ultrarelativistic behaviour of a rotating dust cloud. For this reason, 
we consider the relativistic generalization of the classical Maclaurin disk, which is the flat- 
tened (two-dimensional) limit of the famous Maclaurin spheroids. The corresponding global 
solution to the Einstein equations has recently been found as the solution of a boundary 
value problem 0, first formulated and approximately solved by Bardeen and Wagoner 
13]. In this connection, the appearance of a boundary value problem requires a comment: 
The reason for the complete absence of global solutions describing any uniformly rotating 
perfect fluid ball is that there is no systematic procedure for constructing solutions of the 
non-linear Einstein equations inside the source and matching them to exterior solutions 



along an unknown ('free') surface such that those global solutions are regular everywhere. 
However, in the disk limit, in which the perfect fluid becomes dust^, the body has no longer 
an interior region; it just consists of surface, so to say. Now, the surface conditions (vanishing 
pressure along the surface etc.) can be considered as the boundary conditions on a regular 
stationary and axisymmetric gravitational vacuum field (only the radius remains 'free') and 
this fact makes disk problems accessible to a systematic treatment. Namely, it has been 
shown that 'inverse' methods leading to linear structures (Riemann-Hilbert problems, lin- 
ear integral equations) apply in this case. Following this idea, it was possible to describe 
the gravitational field of the rigidly rotating disk of dust mentioned above in terms of two 
integral equations One of them, the 'small' one, describes the behaviour of the gravi- 
tational field along the axis of symmetry and the physics on the disk (fields, mass-density, 
...). It turned out that its solution can be represented in terms of elliptic functions [H]. 
The other one (the 'big' integral equation) makes use of the solution of the 'small' integral 
equation and describes the gravitational field everywhere. Surprisingly, this equation and its 
corresponding Riemann-Hilbert problem could be solved in terms of hyperelliptic functions 

m 

In the next section we will briefly repeat the mathematical formulation of the boundary 
value problem and its solution. The main intent of this paper, however, is the discussion of 
the physical properties of the disk of dust solution in section 3. Finally, general conclusions 
form the last section. 



2 The boundary value problem and its solution 

To describe a rigidly rotating disk of dust we use cylindrical Weyl-Lewis-Papapetrou coor- 
dinates 

ds^ = e-^^[e^\dp^ + dC) + V] - e^'^idt + ad^f, (2.1) 

which are adapted to stationary axisymmetric problems {U = U{p,(), = cl{p,C)j ^ = 
k{p,())- In these coordinates, the vacuum Einstein equations are equivalent to the Ernst 
equation 

m){f,PP+f<C+-/,p) = fi+fK (2-2) 

for the complex function 

/(p,C)=e^^ + z6 (2.3) 

with 

a,p = pe"^^6,^ , a,^ = -pe~^^b,p (2.4) 

and 

k,p = p[U,l -U,l +\e-'''{b,l -b,l )], k,^ = 2p(f/„ U,^ +^e-^^fe„ b,^ ). (2.5) 

^Strictly speaking, the pressure - to - mass-density ratio tends to zero everywhere. (The maximum 
pressure - at the center - remains finite while the volume mass-density becomes infinite.) The radial 
distribution of the surface mass-density as well as the exterior gravitational field of this object become 
identical with those of a 'genuine' dust disk. 
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As a consequence of the Ernst equation ()2.2[) . the integrabihty conditions 0-,^^ — (^iCpi ^ipC — 
/c,(p are automatically satisfied and the metric functions a and k may be calculated from the 
Ernst potential /. Thus, it is sufficient to consider the Ernst equation alone. 

The metric ()2.1|1 allows an Abelian group of motions G2 with the generators (Killing 
vectors) 

C = Cii < (stationarity), 

7]^ = 6lp, rfrji > (axisymmetry) , (2.6) 

where the Kronecker symbols 51 and 5^ indicate that ^* has only a t-component (,^* = 1) 
whereas r^* points into the azimuthal 99-direction (its trajectories have to form closed circles!). 
By applying ()2.fi|l we get from 1)2.11) the invariant representations 

e'^ = -e.r, a = -e-2%,e (2.7) 

for the "Newtonian" gravitational potential U and the "gravitomagnetic" potential a. 

To formulate the boundary value problem for an infinitesimally thin rigidly rotating disk 
of dust with a coordinate radius po let us start from a spheroid-like rotating perfect fiuid 
configuration (Fig. 1, left) and interpret the rigidly rotating disk of dust as an extremely 
fiattened limiting case of that configuration (Fig. 1, right). 




Fig. 1 — The rigidly rotating disk of dust (right-hand side) as the extremely flattened Umit of a 
rotating perfect fluid body (left-hand side). 

The hydrodynamics of a perfect fiuid with the energy-momentum tensor 

Tik = (e + p)uiUk + PQik, (2.8) 
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e, p, Ui, Qik being energy density, pressure, four-velocity, and metric, respectively, follows 
from the balance equation 

T^';fc=0, (2.9) 
where the semicolon denotes the covariant derivative. Rotational motion of the fluid means 

u' = e-^{C + ^1]'), u'ui = -l, (2.10) 

i.e., the four-velocity is a linear combination of the time-like Killing vector and the 
azimuthal Killing vector r^*. Obviously, 

e^"" = -{C + ^v'm + ^nd- (2.11) 

For rigidly rotating bodies, the angular velocity f2 is a constant, 

Vt = constant. (2-12) 

Then, for an equation of state 

e = e(p), (2.13) 
the pressure p must be a function of V alone, 

p = piV). (2.14) 

This follows from Eqs. ()2.8|) . ()2.9|) . If Eq. ()2.13p is surface-forming, the matching condition 
at the surface (fluid-vacuum interface) requires 

p{Vo) = 0, (2.15) 

i.e. we have 

V = Vo (2.16) 

along the surface of the body. Eq. ()2.11|) inspires us to introduce a corotating frame of 
reference by the transformations 

t' = t, if' = ip- Vtt; (2.17) 

C' = c + ^'n\ V''=V'- (2.18) 
As a consequence of Eq. ()2.18j) . the corotating potential /' constructed from the primed 
Killing vectors with the aid of the Eqs. ()2.7|) and ()2.4|) is again a solution of the Ernst 
equation, 

m')if,pp+f'<c+-/,p) = f'i+f'l (2.19) 

Hence, Eq. ()2.16|) tells us that the real part of the corotating Ernst potential is a constant 
along the surface of the body, 3?/' = exp(2Vo). We may assume the validity of this result 
for the disk limit, too, 

3f^/'lc=o± = e'""" (0<p<Po), (2.20) 

where C = 0^ means top and bottom, respectively. Thus we have found our first boundary 
condition along the disk (above and below). As a consequence of the Einstein equations with 
the perfect fluid source (|2.8p we may conclude in the disk limit that 

(p-ie^^'a'„)„+(p-ie^^'a',c),c = (2.21) 
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holds everywhere including the disk. Applying the procedure known from the transition 
conditions in electrodynamics, one obtains 



C=o+ 



da' 



(2.22) 



C=o- 



at all points of the surface. On the other hand, the reflectional symmetry of the gravitomag- 
netic potential a', 

a'{p,C)=a'{p,-C) (2.23) 



leads to 



da/_ 



C=o+ 



(2.24) 



C=o- 



on the disk, so that da'/d( has to vanish on the disk, 

= (0<p<po). 



(2.25) 



C=o± 



Combining this relation with Eq. ()2.4|) . one has b' = constant on the disk and, after a 
normalization, 

lc=o± = ^'lc=o± =0 (0 < P < Po) (2.26) 

as the other boundary condition along the disk (above and below). Since at infinity (p^+C^ ~* 
oo) the space-time of any isolated source is Minkowskian, we have to take care of the relation 



/I 



(2.27) 



disk: 



infinity: 
/-I 



Po 



Fig. 2 — Boundary value problem [5j. /'is the Ernst potential in the corotating frame of reference 
defined by p' = p, C = C, tp' = — Qt, t' = t {u^' = e~^°6]l). The solution /(p, Q has to be regular 
everywhere outside the disk. 
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The conditions ()2.20|) . ()2.26|) and ()2.27p form our system of boundary conditions for an 
extremely flattened (infinitesimally thin) rigidly rotating disk consisting of a perfect fluid. It 
can be shown that the 2-dimensional mass elements of this disk move on the geodesies of their 
own field. This justifies the designation "disk of dust" , cf. the footnote 1. Fig. 2 illustrates the 
boundary conditions on /. Moreover, the wanted Ernst potential satisfying those conditions 
has to be regular everywhere outside the disk. The axis of symmetry requires particular 
care, especially the elementary flatness condition along the axis (/,p {(,p = 0) = 0). 

The boundary values depend on the two parameters exp(2Vo), cf. fl2.2U|) . and fl, cf. ()2.18|) . 
so that a 2-parameter solution can be expected from the very beginning. Q and exp(2Vo) 
are "source" parameters with a clear physical meaning: Q is the constant angular velocity of 
the mass elements (as measured by an observer at infinity), and exp(2Vo) defines the relative 
redshift 



of a photon emitted from the center of the disk. Obviously, other parameter pairs can be 
chosen. It will turn out that the 'centrifugal' parameter 



is a good measure for the relativistic behaviour of the disk (/i varies between /x = 0: 
Minkowski space and n = Hq = 4.62966...: ultrarelativistic case). We shall present the 
Ernst potential / in terms of yU and Q or fi and po- On the other hand, the solution may also 
be characterized by the total mass M and the ^-component of the angular momentum J, 
which far field quantities can be read off from a multipole expansion. Hence, a connection 



comparable with the parameter relations of "black hole thermodynamics" must hold. 

Let us now briefly outline the solution procedure. We have made use of the so-called 
inverse scattering method of soliton physics, which was first utilized for the axisymmetric 
stationary vacuum Einstein equations by Maison [8 , Belinski & Zakharov [9 , Harrison [Til] . 
Neugebauer [TT], Hauser & Ernst |il2j, Hoenselaers, Kinnersley & Xanthopoulos [12], and 
Aleksejew ^3]. Some of these authors followed the line of Geroch ^15j, Kinnersley [T^ . 
Kinnersley & Chitre J7], and Herlt ^Hl- We have applied a local version ^H] in which the 
Ernst equation (j2.2j) is the integrability condition of the "linear problem" 



Zq = 



(2.28) 



(2.29) 



Q = n{M,J), Vo = Vo{M,J) 



(2.30) 




(2.31) 




(2.32) 



where $(z, ^, A) is a 2 x 2 matrix depending on the spectral parameter 




{K a complex constant) 



(2.33) 
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as well as on the coordinates z = p + iC,, z = p — iC,, whereas C, D and the complex conjugate 
quantities C, D are functions of z, z (p, Q alone. Indeed, from = $,22 and the formulae 

A„=A(A^-l), A„=^(A^-1) (2.34) 

it follows that a certain matrix polynomial in A has to vanish. This yields the set of first 
order differential equations 

C,,= C{D-C)-^{C + D), D,,= D{C-D)-^{D + C) (2.35) 
plus the complex conjugate equations. The system ()2.35|) has the "first integrals" 

Eliminating C and D in ()2.35|) one arrives at the Ernst equation ()2.2|) . Vice versa, if / is a 
solution to the Ernst equation, the matrix $ calculated from fl2.31|l . ()2.32|1 does not depend 
on the path of integration. The idea of the "inverse methods" is to discuss $, for fixed values 
of holomorphic function of A and to calculate C and D from $ afterwards. This 

is an 'inverse' procedure compared with the 'normal' way which consists of the solution of 
differential equations with given coefficients. The term 'scattering' comes from the solution 
technique of the Korteweg-de Vries equation developed by Gardner, Greene, Kruskal & 
Miura [20] whose linear problem has partially the form of the (time-independent) Schrodinger 
equation. In this case, the calculation of the coefficients consists of the construction of the 
Schrodinger potential from the scattering data. 

To construct $ as a function of A we have integrated the linear problem ()2.31|) . ()2.32j) 
along the dashed line in Fig. 1 (right-hand side) and exploited the information of C, D 
and A along the axis of symmetry {f,p= 0, A = ±1), the boundary conditions on the disk 
which simplify C and D' (i.e. we had here to switch to the corotating system) and the simple 
structure of the linear problem at infinity {C = D = 0). In this way, we could pick up enough 
information to construct ^{z, z, A) completely. (The crucial steps were the formulation and 
solution of a matrix Riemann-Hilbert problem in the complex A-plane.) The linear problem 
()2.3H) . ()2.32|) tells us that $ at A = 1 may be normalized in a very simple way, 

^z,z,X = l)=(^ (2.37) 

Hence, once $ is known, the Ernst potential / can be read off from $ (/ = $21 (-2? ^) 1))- 
The result for our problem is [7] 



"""x^dx ""r'x^dx rhx^dx 

+ 



W J W J Wi 

Xl X2 -i 



(2.38) 



where the lower integration limits Xi, X2 are given by 



2 ^ Wf^ ^2 ^ _i±i^ ^^^^ ^ s)^J^^ ^ ^2.39) 
p p 
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whereas the upper hmits Xa, Xb must be calculated from the integral equations 
^"dX ^rdX }hdX ^r^XdX ^f- X dX } hXdX 



W ^ I W J Wi' J W ^ J W J Wi ' ^^'^^^ 

Xi X2 -i Xi X2 -i 

Here we have introduced the abbreviations 

W = WiW2, Wi = ^{X - C/po)' + (p/po)^ W2 = ^lVi^{l + X^ (2.41) 

and 

^ in(Vi + P^(i+xT + Mi + x^)) 

7rzy^l+p2(i + X2)2 

The third integral in ()2.38p as well as the integrals on the right-hand sides in ()2.40|) have to 
be taken along the imaginary axis in the complex X-plane with h and and Wi fixed according 
to 'iRWi < (for p, C outside the disk) and 'Rh = . The task of calculating the upper limits 
Xa, Xb in (2.39) from 

fhdX fhXdX 



is known as Jacobi's famous inversion problem. Gopel |2l] and Rosenhain were able 
to express the hyperelliptic functions Xa{u,v) and Xh{u,v) in terms of (hyperelliptic) theta 
functions. Later on it turned out that even the first two integrals in (j2.38j) can be expressed 
by theta functions in u and v\ A detailed introduction into the related mathematical theory 
which was founded by Riemann and WeierstraB may be found in The repre- 

sentation of the Ernst potential ()2.38|) in terms of theta functions can be found in Stahl's 
book, see [21], page 311, Eq. (5). Here is the result: Defining a theta function {}{x, y; p, q, a) 

00 00 

^{X,y;p,q,a) = ^ ^_i^m+npm^ ^n^ ^2mx+2ny+4mna ^2.44) 

m=— oo n=— oo 

one can reformulate the expressions ()2.38p . ()2.40|) to give 

y _ ^{^qU + aiv - Ci, PqU + piv - C2; p, q, a) ,^ 
{}{aoU + aiv + Ci, PqU + Piv + C2;p,q,a) 

with u and v as in ()2.43|1 and 



} hX^dX 

—i 

The normalization parameters cto, ai, Pq, Pi, 70, 71, the moduli p, q, a of the theta function 
and the quantities Ci, C2 are defined on the two sheets of the hyperelliptic Riemann surface 
related to 



W = fi^iX - Xi)(X - Xi)(X - X2){X - X2){X - tz/po){X + tz/po), (2.47) 
see Figure 3. 





1 1 



where the bracket indicates the characteristic, see |24| . 
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Fig. 3 — Riemann surface with cuts between the branch points Xi and Xi, X2 and X2, —iz/po 
and iz/po- Also shown are the four periods Oj and bi {i = 1,2). (Continuous/dashed hnes belong 
to the upper /lower sheet defined by — > it/iX^ as X — > 00.) 

There are two normalized Abelian differentials of the first kind 

dX XdX 

^dX ^ XdX 
d^. = /5o-^ + A-^ (2.49) 

defined by 

dun = 7ii 5mn {m = 1,2] u = 1, 2) . (2.50) 



Eq. ()2.5U|) consists of four linear algebraic equations and yields the four parameters ao, ai, 
Po, Pi in terms of integrals extending over the closed (deformable) curves oi, 02. It can be 
shown that there is one normalized Abelian differential of the third kind 

dX XdX X^dX 
du^ = 70-^ + 71-^ + /^-^ (2.51) 

with vanishing a— periods, 

j^duj = (j = l,2). (2.52) 
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This equation defines 70, 71 (again via a linear algebraic system). The Riemann matrix 

(^■i)=('2a ('^l.Sl^'^l.a) (2-53) 

(with negative definite real part) is given by 

Bij = jdujj (2.54) 



and defines the moduli p, g, a of the theta function ()2.44|1 . Finally, the quantities Ci, C2 
can be calculated by 



C, = - J dui {i = 1,2), (2.55) 

-iz/po 

where + denotes the upper sheet. Obviously, all the quantities entering the theta functions 
and the exponential function in ()2.45|) can be expressed in terms of well-defined integrals 
and depend on the three parameters p/po, C/Po, P- The corresponding "tables" for ctj. Pi, 
7i, Ci, Bij, u, V, w can easily be calculated by numerical integrations. Fortunately, theta 
series like fl2.44|l converge rapidly. For < /i < po, the solution ()2.45j) is analytic everywhere 
outside the disk - even at the rings —iz/p^ = Xi, X2. 

Figures 4 and 5 give an impression of the Ernst potential for p = 3, as an example. 




Fig. 4 — The real part (e^^) of the Ernst potential for p = 3. 
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Fig. 5 — The imaginary part (b) of the Ernst potential for ji = 2>. 



3 Physical properties 



In dimensionless coordinates, p/po, C/Po, the solution depends on the single parameter p. 
The original parameter Vq entering our boundary value problem can be calculated from p 
via 3?/(p = 0, C = 0"*") (on the axis p = we have 3?/' = 3?/) leading to [U] 

Fo = - - sinh- V p + 5 — i, (3.1) 

2 \^ p[/(p);|p2-4,|p(l + ^)]-|p/' 



/(p) = i r H^ + ^^')d-^ (3.2) 



where is the WeierstraB function defined by 

dt 



X. (3.3) 



Fig. 6 shows exp[2Vo(p)] in the range < p < po = 4.62966184. . . with po being the first 
zero of the denominator in 1)3.11) . This corresponds to > Vq > — oo, where |Vo| -C 1 is the 
Newtonian limit. Note that, according to ()2.29p and ()3.1|) . fipo is a given function of p, and 
we can use either po or VL as the second parameter. 
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Fig. 6 — The corotating Ernst potential /' = e^^° in the disk as a function of the parameter /j,. 

In Fig. 7 the invariant baryonic (proper) surface mass-density cxp can be found as a function 
of p/po for several values of p,. Note that the volume mass-density e entering the energy- 
momentum tensor (Tj^ = eUiUk) may be expressed formally by 



e = a,(p)e^-'=<5(C), 
where S{Q is the usual Dirac delta-distribution. 



(3.4) 



0.20 



0.15 - 



ap/n 



0.10 



0.00 




Fig. 7 — The surface mass-density ap. The normahzed quantity ap/Q is shown as a function of 

P/PO- 
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(Tp can be calculated from 



2% 



(3.5) 



C=o+ 



A more explicit expression has been given in j|6j . From ap we can calculate the total baryonic 
mass Mo, the gravitational mass M, and the total angular momentum J: 



Mn 



Po 



k — U 7 

cTnC pap, 



(3.6) 



po PO 

M = 2j{Tab- -TgabK^dV = 2n J ape"-'' pdp + inQe-'"'' J apc'-^'u^p dp, (3.7) 



PO 

J = - j Tabrfrf'dV = 27re~^° J apc"-^ u'r^ipdp, (3.8) 

E 

where S is the spacelike hypersurface t = constant with the unit future-pointing normal 
vector n". Note that M = exp{Vo)Mo + 2QJ, cf. |2E]- Alternatively, M and J - as the first 
gravitational multipole moments - may be obtained from the asymptotic expansion of the 
Ernst potential, e.g. on the symmetry axis. The dependence of Mq, M and J on p can be 
seen in Fig. 8. 




Fig. 8 — Baryonic mass Mq, gravitational mass M and angular momentum J. The normalized 
quantities JIMq, QM and O^J are shown in dependence on p. 

The relative binding energy (Mq — M)/Mq as well as the characteristic quantity M^/J are 
shown in Fig. 9. 
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1 .0 



0.8 - 



0.6 - 




Fig. 9 — Relative binding energy and M"^ / J . 

For /X — i> /iQ, M and J together with all the other multipole moments approach exactly the 
values of the extreme Kerr solution (^7 = 1/ 2M may be identified with the angular velocity 
of the horizon) In fact, the solution becomes identical with the extreme Kerr solution 
for all values of p and C, except p = C = 0, which represents the horizon of the extreme Kerr 
black hole. Note that, for non-vanishing VL (finite M), po ^ as p ^ po- 



1.4 -_ 

1.2 -_ 

1.0 -_ 

0.3 -_ 
2?lU : 
0.6 -_ 

0.4 -_ 

0.2 -_ 

0.0 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Fig. 10 — Relation between QM and M'^ / J for the classical Maclaurin disk (dashed line), the 
general-relativistic dust disk and the Kerr black-hole 0. 
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Fig. 10 combines the parameter relations between QM and M'^/J for the dust disk as well 
as for the Kerr black-hole. Both branches are connected at the point M^/J = 1. 

Another limit of the space-time for /i — »• /xq is obtained for finite values of p/po and 
C/po, see also The interpretation of the solution for p > po is beyond the scope of 
this paper. We only want to mention here, that there are further zeros pn (n = 1,2, . . .) 
of the denominator in (jH.lj) leading always to the extreme Kerr metric {pi = 38.70908 . . ., 
P2 = 176.92845..., ...). 

A characteristic feature of relativistically rotating bodies are dragging effects due to 
the gravitomagnetic potential. Dragging effects near the rigidly rotating disk of dust have 
been discussed in j^. In particular, the dust disk generates an ergoregion for p > p^ = 
1.68849. . ., see Fig. 11. In this region the Killing vector becomes spacelike. As a con- 
sequence, dif/dt > must hold for any timelike worldline there. Thus, seen from infinity, 
any observer inside the ergoregion is forced to rotate in the same direction as the disk. For 
p ^ Pq the well-known ergosphere of the extreme Kerr black hole appears. 



P = Pe 



p = 1.8 



0( 



p = 2.0 



Of 




fip °= p = 3.0 




P = Po 



Fig. 11 — The ergoregion |31j . 

Note that the shape of the ergoregion for p = 3 may be rediscovered in Fig. 4. 
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4 Conclusions 



In this article, we have considered the rigidly rotating disk of dust as an extremely flattened 
rotating perfect fluid body. Such controllable limiting procedures reducing the dimension of a 
body (here the thickness) are important for the derivation of equations of motion of particles 
from the dynamics of extended bodies. In our case, the motion of the two-dimensional mass 
elements is generally geodesic and independent of the underlying perfect fluid model. In 
that sense, our disk of dust, like the classical Maclaurin disk, represents a "universal" limit 
for any rigidly rotating perfect fluid ball. There is faint hope of an explicit global solution 
for three-dimensional rotating perfect fluid sources with the limit ()2.45|l ! 

Another aspect of our solution is its derivation from a boundary value problem by means 
of the inverse scattering method. It may be expected that this method will prove to be a 
powerful tool for the solution of other boundary value problems for axisymmetric stationary 
gravitational vacuum fields in Einstein's theory. It could also improve the insight into the 
structure of the exterior fields of rotating bodies. A first step in this direction consisted in 
the treatment of the reflectional symmetry of the gravitational field with the aid of the linear 
problem (imT|l . (jZ!?^ in 

Finally, as a 'practical' application, the solution could be used as a testbed for numerical 
codes describing rotating star models in general relativity, as, e.g., neutron stars. 

There is also a formal aspect to be mentioned. Generalizing the expressions fj2.38|) . 
()2.40|) it was possible to construct a solution class in terms of hyperelliptic theta functions 
containing an arbitrary potential function and - depending on the genus - an arbitrary 
number of constants ^28^. The relation of this class to the finite-gap class of solutions [221 
requires a subtle discussion of certain limiting procedures. We want to emphasize that it is 
not possible to know a priori into which class of (known or unknown) solutions of the Ernst 
equation a boundary value problem might fall. 
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